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We present the constraints on neutrino oscillation parameters Am Q and #0 governing the solar neutrino oscilla- 
tions from two generation analysis of solar and KamLAND data. We include the latest 766.3 ton year KamLAND 
data in our analysis. We also present the allowed values of parameters AmJ im and sin 2 9 a tm from two generation 
oscillation analysis of SuperKamiokande atmospheric and K2K data. For both cases we discuss the precision 
achieved in the present set of experiments and also how the precision can be improved in future. We also obtain 
the bounds on #13 from three generation analysis of global oscillation data. We emphasise on the roles played by 
different data sets in constraining the allowed parameter ranges. 



1. Introduction 

Compelling evidence in favour of neutrino os- 
cillation first came from atmospheric neutrino 
flux measurements at SuperKamiokande (SK) pQ. 
The solar neutrino shortfall at SAGE, Gallcx, 
GNO 2 , Kamiokande and SK |3] were sugges- 
tive of neutrino oscillation. This was firmly estab- 
lished by SNO @E] and KamLAND (KL) [oT7] . 
The later provided the first evidence of oscilla- 
tion of reactor neutrinos. The K2K experiment 
demonstrated oscillations of accelerator neutrinos 
8 . Earlier the LSND 9 experiment reported 
positive evidence of (y^) oscillation using low 
energy accelerators though this was not confirmed 
by the similar experiment KARMEN. The MINI- 
BOONE experiment will provide an independent 
check |10j . 

Aim of global oscillation analysis is to analyze 
the experimental data using a suitable statistical 
procedure and extract the information on neu- 
trino oscillation parameters Am 2 (mass squared 
difference in vacuum) and 9 (mixing angle in vac- 
uum). Various statistical methods have been 
adopted starting from the standard covariance 



method, pull method |ll|.Freauentist method 22 
and Bayesian Analysis ^3] ■ 

2. Solar Neutrino Oscillation parameters 

In this section we discuss the constraints on 
solar neutrino oscillation parameters Am 



Amli, 9 



analysis. 



9\2 from two flavour v e 
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Allowed area from global Solar Data : In 

our analysis of global solar data we include the 
total rates from the radiochemical experiments CI 
and Ga (Gallex, SAGE and GNO combined) 
and the 1496 day SK Zenith angle spectrum data 
H3, the SNO spectrum data from the pure D%0 
phase [U, as wen as the CC, NC and ES rates 
from the salt phase [5]. We use BP04 fluxes [T4*j 
and keep the S B flux normalisation free. Further 
details of our analysis can be found in [15116117] . 
In the left panel of Figure ^ we show the allowed 
region in the parameter space including the SNO 
spectrum data from the pure D%0 phase while 
the right panel shows the allowed regions after 
including the data from the salt phase. The best- 
fit values of parameters obtained are 
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• Am| x =6.1xlO- 5 eV 2 , sin 2 12 = 0.3, / B =0.89. 

•99% C.L. range (before salt) 

Am^ =(3.1-25.7) xlCT 5 eV 2 , sin 2 12 =O.21 -0.44 

•99% C.L. range (after salt) 

Aml x =(3.2-14.8) xlO- 5 eV 2 , sin 2 6>i 2 =0.22 -0.37. 

The allowed region corresponds to Large-Mixing- 

Angle (LMA) MSW solution. The impact of 



Solar(BP04) [pre-salt] 


Solar(BP04) [post-salt] 


0.6 o.5 




& 


0.6 0.5 


0.5 


0- 




0.45 


0.45 




0.2 0.3 0.4 


0.2 0.3 0.4 



0.1 0.2 0.3 0.4 0.1 0.2 0.3 0.4 0.5 

sin 9 12 



Figure 1. The 90%, 95%, 99% and 99.73% C.L. 
allowed regions in the Am^ —sin 2 9q plane from 
global x 2 -analysis of solar neutrino data. Also 
shown are the iso-CC/NC contours. 



the salt data is to tighten the upper limit on 
Am| x and sin 2 The iso-CC/NC curves in 
Fig. ^ shows that for higher values of sin 2 9\2 
and Am^ the CC/NC value is higher. Since 
the CC/NC ratio for the salt phase data is 0.31 
as opposed to 0.35 in the pure D 2 phase the al- 
lowed regions shift to smaller sin 2 (?i 2 and Am 21 
values following the iso-CC/NC contours. 

Allowed area from KamLAND spectra: 

We next study the impact of the new results 
from the KL experiment in constraining the os- 
cillation parameters. The survival probability for 



162 Ton-yr KL data 766.3 Ton-yr KL data 




tan 9 12 



°~ Figure 2. Same as in figure 1 but in in the 
Amji — tan 2 9i2 plane from global x 2-ana ly s i s 
of KL spectral data. 



KL is : 

5 kt 9 9 / 1.27Am 2 1 L\ , , 

P« L = 1 - sin 2 29 12 sin 2 f (1) 

neglecting the small matter effect for lower values 
of Am|i. Assuming CPT conservation i.e same 
oscillation parameters for neutrinos and antineu- 
trinos KL is sensitive to Am^i in the LMA re- 
gion. 

In the right panel of fig [21 we show the al- 
lowed regions obtained after including the latest 
KL data in our analysis while the left panel shows 
the allowed regions without this data. The figure 
shows that higher Am 2 regions, which correspond 
to an averaged survival probability and hence flat 
spectra, reduce in size due to the increased dis- 
tortion in the new KL spectrum data jZj. For 
details of our KL analysis we refer to |18ll9j . 

Allowed area from Solar +KamL AND: In 

left (right) panel of|3|we present the allowed re- 
gions from analysis of solar and 162 (766.3) ton 
year (Ty) KL data. In the left panel there are 
two allowed zones - the low-LMA zone which is 
more favoured and the high-LMA zone admit- 
ted only at 3cr. In the right panel there is no 
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Figure 3. Same as in figure 1 but from global 
X 2 -analysis of solar and KL data. 



Figure 4. A\ 2 vs Am^ (sin 2 6*12) marginalised 
over the other parameters. 



Am 2 ! 
Ami 



allowed high-LMA region at 3a since the latest 
KL data disfavour higher Am| x values. The 
allowed range of Amfi in low-LMA region also 
narrows down considerably. The best-fit points 
from combined analysis of solar with 162 Ty and 
766.3 Ty KL data are found respectively to be 
7.2 xl(T 5 eV 2 , sin 2 0i 2 = 0.29, 
, - 8.3 xlCT 5 eV 2 , sin 2 9 12 = 0.27. 
The 9 > 7r/4 (Dark-Side) solutions allowed by 
KL spectrum data get disallowed by the so- 
lar data as the matter effect in sun break the 
6 — > 7r/2 — 9 symmetry (cf. fig. |SJ). 

In fig H we plot Ax 2 = X 2 - x 2 min vs Am 2 21 
(right panel) and vs. sin 2 9 12 (left panel). 
This figure shows that with increased statistics 
KL data constrains Am|x more and more sharply. 
However sin 2 #12 is not constrained as much. 
Maximal mixing is found to be more disfavoured 
(> 5a) by the new KL data because of increased 
precision. 

The figure [S] shows the probabilities vs energy 
for KL (at an average distance of 180 km) and so- 
lar neutrinos for Am| x =7x 10~ 5 eV 2 and sin 2 12 
=0.3. For these parameters the high energy 8 B 
solar neutrinos undergo adiabatic MSW transi- 
tion with P ee ps sin 2 9. For low energy solar neu- 
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Figure 5. The probabilites for solar and KL as a 
function of energy 
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trinos P ee w 1 — 0.5 sin 2 20. Thus with decrease 
in energy the survival probability increases for 
9 < 7r/4 causing an upturn in the observed en- 
ergy spectrum. Whereas for the solar probabili- 
ties the Am|i dependence is completely averaged 
out in KL the probability exhibits a L/E depen- 
dence which gives it an unprecedented sensitivity 
to Am 2 ! . 

The imprint of MSW effect |U] in Sun in 
the solar and KL data can be quantified by 
parametrising the matter induced potential in 
the sun, Vmsw — \/2GFn e , as ausw-V |21j. 
ctMSW — corresponds to no MSW effect while 
(xmsw = f corresponds to full MSW effect. 
Analysis performed in |21| with solar and f62 
Ty KL data shows that cxmsw = case is re- 
jected at 5.6er w.r.t the minimum which comes 
at o-msw =1 i- e for standard MSW effect. How- 
ever the allowed range of umsw was found to 
be large. Increased statistics from KL can put 
stronger constraints on otMSW and can provide 
more precise test of MSW and "new" physics 
beyond MSW. 
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Figure 6. The allowed regions in Am^ -sin 2 612 
plane, obtained by removing the data from one so- 
lar neutrino experiment from the global fit. These 
figures are drawn on a linear scale for Am^ . 



Impact of each solar experiment: In Fig. 
El we present the allowed regions obtained by 
taking out the data from one solar neutrino ex- 
periment from the global data set. Taking out the 
SNO data allows smaller values of fs, and hence 
larger values of of sin 2 9\2 (P ee ~ /s sin 2 #12). 
With f 62 Ty KL data excluding SNO allowed 
maximal mixing, Dark side and higher Am 2 so- 
lutions. With the spectral distortion observed in 
the 766.3 Ty KL data these regions are ruled out 
with increased confidence even without SNO. 

On the precision of v Q oscillation parameters: 



set of solar experiments. The reason for this is the 
average energy and distance in KL corresponds 
to a Survival Probability MAXimum (SPMAX) 
i.e sin 2 (Am| 1 L/AE) » 0. This means that the 
coefficient of the sin 2 29±2 term in P^ L is rela- 
tively small, weakening the sensitivity of KL to 
#12. The precision in Qyi can be improved by re- 
ducing the baseline length such that one gets a 
minimum of the v e survival probability (SPMIN) 
where sin (Ato^i L/4E) = 1. This corresponds 
to a distance L = 1.2A(E/MeV)(eV 2 / Amf 2 ) m. 
For the low-LMA solution region and the average 
energy of neutrinos in KL experiment, this corre- 
sponds to a distance of approximately (50 - 70) 
km |23| . For an experiment with a 70 km base- 
line and 24.3 GW reactor power sin 2 612 can be 
determined with ~ f 0% error at 99% C.L. with a 
3 kTy statistics [23j . 

Improved sensitivity in the measurement of 612 
is possible from LowNu experiments 22 . The pp 
flux is known with 1% accuracy from Standard 
Solar Models as compared to the uncertainties of 
S B ~ 20% and 7 Be ~ 10%. At the low energies 
relevant for pp neutrinos P ee ~ 1 ~ 5 sin 2 29 re- 



in Table [1J we present the 3<t allowed ranges of 
Ami^ and sin 2 #12, obtained using different data 
sets. We also show the uncertainty in the value 
of the parameters through a quantity "spread" 
which we define in general as 

spread - Am2 ( sin2 8) max ~ ATO 2 (sin 2 9) mm 
Am 2 (sin 2 9) mra + Am 2 (sin 2 9) mm 

Table I illustrates the remarkable sensitivity of 
KL in reducing the uncertainty in Ato^i- But #12 
is not constrained much better than the current 
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Table 1 



3d allowed ranges and % spread of Ato^i an d sin 2 612 obtained from 2 parameter plots. 



Data set 


(3<r)Range of 


(3cr ) spread in 


(3a) Range of 


(3cr) spread in 


used 


Am 2 21 eV 2 


Amfi 


sin 2 6»i2 


sin 2 612 


only sol 


3.0 - 17.0 


70% 


0.21 - 0.39 


30% 


sol+162 Ty KL 


4.9 - 10.7 


37% 


0.21 - 0.39 


30% 


sol+ 766.3 Ty KL 


7.2 - 9.5 


14% 


0.21-0.37 


27% 



> 8 



S + K 3 yr + [ 7 Be] v e ± 5% 
+ [p-p] v . e + 3% 
+ [p-p]„±l% 



0.3 




0.4 



0.5 



tan 2 9, 



0.6 



Figure 7. Allowed area including the rates from 
a generic pp neutrino experiment with different 
errors. Figure taken from |24| 




suiting in a heightened sensitivity to 9\2- In|J|we 
show the allowed region [23] in Am 2 — tan 2 612 
plane including fictitious results from a generic 
pp scattering experiment. The figure shows that 
precision in 612 increases with reduced error in pp 
flux measurment. 

The recently proposed Gd loaded SK detector 
can measure Ato^i with ~ 1% and sin 2 812 with 
~ 15% uncertainty at 99% C.L. after 5 years of 
data taking [2*5], 

3. Atmospheric Neutrino Oscillation Pa- 
rameters 

In this section we present the constraints 
on atmospheric neutrino oscillation parameters 
Am§ 2 = Am 2 tm and sin 2 O23 = sin 2 6 atm from 



Figure 8. The allowed regions from SK (line con- 
tours) and SK+K2K (shaded contours) analysis 
from a two flavour analysis. Also shown is the 
A% 2 vs Am 2 tm and sin 2 9 atm marginalised w.r.t 
the undisplayed parameter|28| 



two flavour oscillation analysis SK and SK+K2K 
data. 

Allowed parameters from SK+K2K analysis: 

In 2003 the SK collaboration reanalyzed their 
data and reported a downward shift of the al- 
lowed Am 2 tm from a two-generation analysis 
due to incorporation of three-dimensional Honda 
atmospheric fluxes and revised cross-section and 
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efficiencies. The best-fit value of Am 2 tm from SK 
zenith angle data reported in this conference is 
2.1xlO~ 3 eV 2 |2S]- Analysis of SuperKamiokande 
data by two groups confirm this downward shift 

of AmL Baza. 

The K2K experiment which is sensitive to at- 
mospheric neutrino oscillation parameters gives 
the best-fit as Am 2 tm = 2.8 xlCT 3 eV 2 
sin 2 20 atm = 1.0. 

In fig. [HJwe show the allowed regions in Am 2 tm 
- sin 2 Oatm plane from analysis of SK data as well 
as from combined analysis of SK+K2K data [2Hj • 
The best-fit from the combined analysis corre- 
sponds to 

• Am 2 atm = 2.3 xKT 3 eV 2 , sin 2 8 atm = 0.5. 
Higher Am 2 values are seen to be constrained by 
K2K but Oatm is not constrained any better. For 
a two generation analysis the relevant probabil- 
ity for both SK and K2K are vacuum oscillation 
probability and the allowed regions are symmetric 
about 6*23 and ir/2 — 823- 

Recently SuperKamiokande has reported the 
observation of the first oscillation minima in 
the L/E distribution (2Hj ■ The observation of 
this first dip disfavours the neutrino decay and 
de-coherence solutions. The best-fit value of 
Am 2 tm from this analysis is 2.4 xl0~ 3 eV 2 with 
sin 2 29 atm — 1-0. The allowed ranges of parame- 
ters from the two different analyses are: 
•90% C.L. range (SK L/E) 
Am 2 tm = 1.9 - 3.0 xl0~ 3 eV 2 , sin 2 29 atm > 0.9 

• 90% C.L. range (SK Zenith) 

Am 2 tm = 1.3 - 3.0 xl0~ 3 eV 2 , sin 2 20 atm > 0.9 

Precision of atmospheric neutrino oscillat- 
ion parameters: The spread in Am 2 tm is 39% 
from SK zenith data whereas it is 22% from SK 
L/E data. Thus improved precision in Am 2 tm 
is obtained with L/E data. Range of sin 2 8 a t m 
however remains unchanged: (5(sin 2 2^23) ~ 5%. 
However <5(sin 2 6>23) ~ 32% because sin 2 #23 pre- 
cision is worse than sin 2 2023 precision near max- 
imal mixing. Increased statistics in L/E data 
can improve Am| 2 uncertainty to ~ 10% at 90% 
C.L. in 20 years of operation of SK [301 • This 
is sensitive to the true Am§ 2 chosen and the 
above precision is achieved for Am 2 2 =2.b xlO -3 



eV 2 . But precision in sin 2 #23 does not improve 
with increased statistics jSO]- Large Magnetized 
Iron calorimeters for atmospheric neutrinos (e.g 
MONOLITH,INO) [5T152] can improve the pre- 
cision in both atmospheric parameters. 

4. Three Flavour Oscillation 

The oscillation parameters for a three flavour 
analysis are two mass squared differences Am^ — 



Am 2 ! 



^ m CHOOZ 



Ami 



Am 
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and three mixing angles which are usually de- 
noted as #12, #23 & nd #13 with the mixing matrix 
U parameterized in the standard MNS form. In 
the limit Am\ x « AmJ 2 as is indicated by the 
two generation analysis of atmospheric and so- 
lar neutrino data. In this limit the atmospheric 
probabilites depend on Am 2 2 , 6*13, 623 and so- 
lar neutrino probabilites depend on Am^, 9\%, 
613. The CP violation phases can be neglected 
in this limit. Clearly for 6*13 = solar and at- 
mospheric neutrinos decouple. The probability 
for the CHOOZ reactor experiment depends on 
Amjj and #13. 

Bound on sin 2 9 13 : CHOOZ/PaloVarde ex- 
periment has put stringent bounds on 813 from 
non-observation of i7 e disappearance [HSj- In fig 
El we show the allowed area from CHOOZ experi- 
ment in Ato 2 ^ — tan 2 #13 plane. The area shaded 
by vertical blue (slanting red) lines correspond to 
the allowed range of Am^ from two generation 
analysis of SK L/E (zenith) + K2K atmospheric 
neutrino data |34I35| . The figure shows that the 
813 bound from CHOOZ depends sensitively on 



Ato^. Stronger bounds are obtained for higher 
Am 2 .!. 

In figure ITUI we show the limits on 813 obtained 
from A\ 2 vs #13 plots. We marginalise over 
all other parameters. We let Am| x vary freely 
within the 3cr range obtained using the one pa- 
rameter Ax 2 vs Am]; fit as given in the latest 
SK analysis |2BJ. The bounds that we get are 

• sin 2 0i3 < 0.071 (CHOOZ+atm) 

• sin 2 0i3 < 0.063 (sol+KL (162)+CHOOZ+atm) 2 

• sin 2 0i3 < 0.05 (sol+KL (766.3)+CHOOZ+atm) 



2 Using the rage of Ami from the SK+K2K analysis of 
1521 bound comes as sin 2 6*13 < 0.077 El 



7 




95% CL 

99% CL 

| " I jm CL : 

1(T 10" 2 10"' 10° 10 1 10 2 

tan 2 9, 3 



Figure 9. Allowed area in Am^ — tan 2 #13 plane 
from CHOOZ data. The shaded regions corre- 
spond to the allowed range of Am§ x using SK 
zenith angle (L/E) and K2K data |34I35| . 
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A non-zero #13 prefers higher Am^ values which 
is disfavoured by the solar+KL data. Because 
of this the combined \ 2 increases making the 
bound on #13 stronger than compared to only 
CHOOZ+atmospheric. Since the recent KL data 
disfavours high Am| x values to a greater extent 
the bound on #13 further improves by including 
it. 

In ^3 we show the effect of non-zero sin 2 #13 
on the allowed area in solar neutrino parameter 
space. The presence of a small non-zero #13 can 
improve the fit in the regions of the parameter 
space with higher values of Am|i Eh i- e -i m 
the high-LMA zone. But due to increased dis- 
cord between the new KL data and high Ato^i 
values the high-LMA region now gets excluded at 
more than 3<r even in the presence of a third gen- 
eration in the mixing, indicating the robustness of 
the low-LMA solution. The allowed atmospheric 
regions in Am^- sin 2 2#23 plane also remains sta- 
ble in presence of non-zero #13 |2Hj ■ Including the 
ATO21 terms in the analysis one can give bound on 
the hierarchy parameter a — Am^/ Am 2 ^ associ- 
ated with subleading oscillations for atmospheric 
neutrinos as well as for long baseline studies. The 
best-fit for this comes as a = 0.03 



5. Conclusions 

A series of seminal experiments performed over 
a period of nearly four decades and more and 
more refined analysis of the global experimen- 
tal data has not only confirmed the existence 
of neutrino mass and mixing, but also narrowed 
down the allowed parameter ranges considerably 
heralding the precision era in neutrino physics. 

The prime future goals for neutrino oscillation 
studies are- further improvement of precision of 
the parameters, to obtain some measure of small- 
ness of #13 , to determine the sign of Am| 2 and 
to search for CP violation in the lepton sector. 



Figure 10. Bounds on the mixing angle #13 from 
data combination. Also shown is the 3cr A% 2 limit 
for 1 parameter. 
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Figure 11. Same as in figure 1 but from a three 
flavour analysis of the global solar and reactor 
neutrino data and plotted on a linear scale for 
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